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In this paper we present a preliminary optimal orbit analysis 
for the Laser Interferometer Space Antenna (LISA). LISA is a 
NASA/ESA mission to study gravitational waves and test pre- 
dictions of general relativity. The nominal formation consists of 
three spa.cecraft in heliocentric orbits at 1 AU and trailing the 
Eartli by twenty degrees. This configuration was chosen as a 
trade off to reduce the noise sources that wdll affect the instru- 
ment and to reduce the fuel to achieve the final orbit. We present 
equations for the nominal orbit design^ and discuss seveial dif- 
ferent measures of performance for the LISA formation. All of 
the measures directly relate the formation dynamics to science 
performance. Also, constraints on the formation dynamics due 
to spacecraft and instrument limitations are discussed. Using 
the nominal solution as an initial guess, the formation is op- 
timized using Secpiential Quadratic Programming to maximize 
the performance while satisfying a set of nonlinear constraints. 
Results are presented for each of the performa.nce measures. 


INTRODUCTION AND MISSION OVERVIEW 


LISA is a NASA/ESA mission to detect and study gravitational waves from massive 
bla.ck hole systems and galactic binaries. Gravitational waves are ripples m space-time 
caused by massive objects. They are a prediction of general relativity but are yet to be 
directly detected. The detection and understanding of gravitational waves can piovK e 
breakthroughs and refinements in current relativistic theory. 

EISA’s primary objective is to detect gravitational waves in the frequency range 
of to 1 Hz This frequency range is chosen to complement ground observations 

for frequencies more than 100 Hz. The primary reason for using a space-based system 

is to eliminate gravitational noise sources. Ground based observatoiies aie imi e o 

high frequency measurements because of the extreme sensitivity to gravity gradient and 
seismic noise. 


In this paper we begin by discussing the LISA mission primarily from a ni.ssion 
design perspective. Secondly, the science objectives and their implications on the orbit 
design are discussed. We consider several possible measures of performance. Next, we 
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present the nominal orbit and discuss its performance under two-body and real-woild 
orbit dynamics. Using the nominal orbit as an initial guess, a Sequential Quadratic 
Pro<^ramining (SQP) / Projected Gradient approach is used to optimize the formation 
creometry while simultaneously satisfying design constraints. Results are pi^sented for 
Lch of the performance measures considered. Finally, we disc-uss some conclusions and 
point out where more work is necessary to expand the set of acceptable solutions. 


MISSION OVERVIEW 


In this section we giye an oyerview of the LISA mission. only present 


brief 

summary of topics that are interesting from an orbit perspectiye. Comsiderable concep 
deyelopment work has been performed in most areas of the mission, and for moi e detai e 
information on areas not directly related to the orbit design, we refer the reader to the 
LISA Final Technical Report.^ 


The nominal LISA formation consists of three spacecraft in carefully designed orbits. 
The relatiye positions of the spacecraft maintain a nearly equilateral triangle configura- 
tion as seen in Figure 1. The minimum science can be performed with two of the legs je 
the third leg can proyide additional information as well as offer redundancy in the case 
that one leg is lost. The primary measurement is the rehatiye change m length of two of 
the legs of formation. This measurement can giye information about both the dnection 
and polarization amplitude of the grayitational waye. 


Optical Assembly 



Figure 1: Equilateral Formation Configuration 


Each of the three spacecraft has two optical assemblies which can be articulated 
so that each assembly points at the appropriate spacecraft as the formation eyo yes^ 
The variation in the angle of articulation is on the order of a few degrees. Each optical 
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assembly can transmit and rec;eive a laser signal. One significant difference between LIS. 
and a conventional interferometer is that the transmitted signal is not simply leflected 
back. Because the leg lengths are 5 million km there is significant loss of signal strength 
over the traversal. To ensure adequate signal strength, the laser on the receiving end 
is phase locked to the incoming signal, and tlie phase, locked signal is then transmitted 

back. 

The leg length measurement is referenced to a free-floating proof mass inter nal to 
the spacecraft. The spacecraffs Disturbance Reduction System (DRS) is designed to 
ensure that no non-gravitational disturbances affect the proof mass trajecton in the 
measurement direction. The DRS is also responsible for ensuring that the spacecraft 
point in the appropriate direction. Self gravity between the spacecraft and the proof mass 
can cause undesirable complications in the DRS design and unacceptable measurement 
noise and is a topic of current research. 

The target launch date for LISA is 2010 - 2012. Each spacecraft will follow an inde- 
pendent, optimal low-thrust trajectory to acliieve the mission orbit. Once the spacecraft 
arrive in the mission orbits, they separate from the propulsion modules to reduce any dis- 
turbance that the module might introduce into the measurement. The nominal mission 
life is two years with a possibility for an extended mission of five to ten years. 

ORBIT REQUIREMENTS AND PERFORMANCE 

In this section we discuss several different sets of orbit performance metrics that 
are of interest for LISA. The orbit requirements are driven by the signal processing 
approach. We discuss five performance measures here and present results for each in a 
later section. Part of the challenge of the LISA orbit optimization is finding an optimal 
solution within the mission design constraints. Wn conclude this section with a discussion 
of some nonlinear constraints that each solution must satisfy for a feasible orbit. 

The current approaches for LISA signal processing can be broken down into two 
categories. Here we present a high-level view of the approaches. For a detailed discussion 
see Tinto.^ The first category uses data from two of the legs in the formation. Two legs 
of the interferometer can provide the minimum science. In this case the third leg provides 
redundancy and additional measurement information. The second category involves using 
measurements from all three legs in the signal processing. For either category of signal 
processing, Doppler shifts due to relative velocities above about 15 m/s can degrade 
performance. So, it is necessary to ensure that the maximum rate-of-c:hange of any leg 
of the formation is less that 15 rn/s. 

Before discussing the individual performance measures, we define a few useful vari- 
ables. The position and velocity vectors of the spacecraft are denoted r, and f 
rsspGctivclv- ThG iGg forniGcl by sp3.cc(‘r8.ft i cind j is denotGd 

L,j=r,-rj (1) 
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a.nd the leg length is simply gi\'en by 


= II r,- - r,l| 

The rat,e-of-change of the leg formed by the and is denoted 

dL,j _ dr, drj 

dt dt dt 

and the magnitude of the rate of change is simply 


of j 


dr,; dvj 

dt 


dt dt 


( 3 ) 


( 4 ) 


Two of the five performance measures considered are category one. The first measuie 
considered, labeled cfj, is based on the rate-of-change of two legs of the formation. By 
minimizing the difference in rate-of-chaiige between two of the legs in the formation we 
may be able to improve the science performance. The first cost function has the foim 



Cl 




dL23 

dt 


dt 


( 5 ) 


where r, is the velocity vector of the zt.h spacecraft, j, and k are indices for the remain- 
ing two spacecraft, and Cj is a constant for normalization. By minimizing c/i , we will 
nfinimize the difference in relative velocities between two of the legs of the formation. 


The second cost function considered is based on the difference m lengths of two of the 
legs of the formation. There are two motivations for a cost function of this form. One 
reason is that equal leg lengths may provide better science. The second reason is that if 
we achieve equal leg lengths over the entire mission, we achieve equal leg rates as well. 
This is beca.use if the leg length functions are the same then so are their derivatives. Ibis 
cost function is labeled c/2 and is expressed as 

c/2 = - ^{L^2-Lnfdt (6) 

C2 Jo 


The remaining three cost functions are category two, and are concerned with opti- 
mizing a metric based on all three of the legs of the formation. The third cost function 
considered is the average rate of change of the legs of the formation. The function is 
labeled c /s and can be expressed as 


= - r 

C3 JO 


dL 


12 


dt 




dL23 


dt 


+ 


dL 


13 


dt 


dt 


It is also of interest to determine if we can find an orbit solution that yields all 
leg lengths of the formation equal in time. This solution, assuming that the design 
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constraints are satisfied, would minimize difficulties in the signal piocessing scheme as 
well as mitigate redundancy difficulties if one the legs of the formation is lost. The fourth 
cost function, c/j. is expressed as 


.r = i ((L,o - Loa)- + {Lu - + (/^i2 - Lnf) dl 

ca Jo ^ 


( 8 ) 


where L19 is the 


> leno-th of the leg formed by spacec:raft one and spacecraft two etc:. 


Finally another approach to ensure that the three legs are ahvays equal in length is 
to require that the interior angles of the formation are always 60 degrees. This would 
give a ‘Treathing-^ equilateral triangle. The final cost function is labeled c/5 and can be 

expressed as. 


cf. = 1 f ' ({dusit) - 7t/ 3)2 + (0231 (i) - ^/3)' + (0312(0 ~ 7 t/3)") dt (9) 

C5 Jo ^ 

where Ons is the angle subtended by the vectors from spacecraft 1 to spacecraft 2 and 
from spacecraft 2 to spacecraft 3. 

The challenge of the orbit optimization is to find solutions that minimize a particu- 
lar cost function and also satisfy a set of nonlinear constraints imposed by instrument 
and communications limitations. As mentioned previously, the largest relative velocity 
between any trvo spacecraft must be less than or equal to 15 m/s. The absolute eg 
lengths are also a concern. The length of the legs determines the frequency range of the 
instrument. The nominal leg length is chosen as 5 million km to ensure the instrument 
is sensitive to the desired frequencies of gravitational waves. We constrain the solutions 
so that the leg lengths he within 2% of the nominal leg length. Due to articulation limi- 
tations on the spacecraft, we constrain the interior angles of the formation to he between 
58.5 degrees and 61.5 degrees. Finally, to meet communications and orbit determination 
requirements, the formation must be no more than twenty degrees away from the Earth, 
as measured from a vertex at the Sun. 


Before moving on to the developement of the nominal formation, it is impoitant to 
note that there are other important measures of performance that are not considered 
here. Other possible measures of performance include the amount of fuel required to 
reconfigure the formation from one optimal configuration to another m the event that 
one of the legs is lost. How much fuel is required to reconfigure from one optimal scenario 
to another is a very important question because the onboard thrusters are expected to 
provide a maximum of about 80/rN of thrust. Hence, reconfiguring the formation could 
take on the order of months, or for some configurations be completely infeasible. 
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THE NOMINAL ORBIT 


In this section we present the nominal orbit to achieve a nearly equilateral configu- 
ration using basic two-body orbit dynamics. The approach used here is based on work 
from Folkner^ et. ai. The nominal orbit does not necessarily satisfy all of the oibit le- 
quiremcnts. However the nominal orbit does provide a good initial guess for a numenc:al 
optimization approach. We discuss the two-body evolution of the nominal orbit as well 
as the evolution of the formation in the presence of real-world perturbations. 

A set of orbital elements which provides a nearly eciuilateral formation is given below. 
Note that the semimajor axis a, the eccentricity, e, the inclination, i, and the argument 
of periapsis. m. are the same for all orbits in the foimation. 


a — a£ 

d 

2 \Zia.E 
d 

uj = 7t/ 2 or 37 t/2 


( 10 ) 

( 11 ) 

( 12 ) 

(13) 

(14) 


where as is the semimajor axis of the Earth’s orbit about the Sun, and d is the nominal 
leg length. The longitude of the ascending node, fl, and the mean anomaly, M are 
different for each orbit and are given by, 

for S/Cl - (a A/) (15) 

for S/C 2 = (Q + 27t/3, A/ - 27t/3) (16) 

for S/C 3 = (L> - 27t/3, M + 27t/3) (17) 


Before investigating the evolution of the nominal orbit, it is useful to define some 
parameters that characterize some important aspects of the motion. The average leg 
lengtli over the mission life is denoted L and is given by 

L= r {Ln + L2,)dt (18) 

2tf Jq 

for a category one cost function and is given by 

L — — [ {Li2 -f L 23 + Ly.^) dt (19) 

Stf Jo 

for a category two cost function. The maximum difference between legs of interest, over 
the mission life, is defined as Ar. Note that this quant it}^ is defined slightly diffeientlv 
depending on whether the cost function is category one or two. For a cost function aimed 
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at optimizing a parameter based on two legs (category one), only the tw-o legs of inteiest 
are used in determining Ar and we have 

Ar = max (/’.i 2 (/) - f'2.3(0) 


For a category two performance metric it is more useful to define Ar as 

Ar = max ( [(fi] 2 (^) ~ ^■^23(0) (^-^i2(0 ~ (^2s(0 ~ ^"i3(f))] ) 


( 21 ) 


Another useful design parameter is the maximum percent difference in leg length. 7oAr 


and is given by 


%Ar = lOOAr/L 


( 22 ) 


The maximum difference in leg rate of change is denoted Ar and for a category one 
metric is given by 

Ar = max 

For a category two metric we have Finally, the maximum leg rate of change is denoted i 
and for a category one metric; we have 



dhi2{i) 


dL23(^) 



dt 


dt 

1 


(23) 


r = max 


dLn^jt) 

(it 


dli 23 {t) 

dt 


(24) 


and for a category two metric we have 


f = max 


dL^t) 


dlj2‘s{i) 

dt 


di 


dLisjt) 

dt 


(25) 


Figure 2 shows the evolution of several relevant parameters of the nominal solution 
under the gravitational influence of a spherical central body with no perturbations. The 
upper left hand plot shows the lengths of the legs of the formation as functions of time. 
The leg lengths oscillate sinusoidally centered at 5 million km with an amplitude of 
about 100 thousand km. In the top right plot we see the rate of change of the leg lengths 
over one year. The leg rates oscillate about zero m/s with a period of one year and an 
amplitude of about 4 m/s. In the bottom left hand plot of Figure 2, we see the evolution 
of the interior angles of the formation. The angles oscillate about 60 with an amplitude 
of about 0.45 degrees. 

In table 1 we see the relevant statistics for the nominal orbit. The average arm length 
is about 5 million km as desired. The maximum difference between any two arm lengths 
is 39370 km and the maximum percent difference is 0.784%. The maximum difference m 
leg rates is 7.948 m/s and the maximum leg rate of change is 4.010 m/s. 

Before moving on to the numerical optimization, we look at the affects of pertur- 
bations on the evolution of the nominal orbit. Figure 3 shows the evolution of some 
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^ Range vs. Time 




Range Rate vs. Time 



Figure 2: Nominal Unperturbed Evohttiou 


Table 1: Statistics for the Unperturbed Orbit 


Variable Name 

2 Leg Definition 

L 

4.999e^ km 

Ar 

39370 km 

%Ar 

0.7874 % 

At 

7.948 m/s 

r 

4.010 m/s 


relevant parameters when perturbations arc included in the orbit force model for a five 
year propagation. All of the planets and the Earth’s moon are included in the force 
model. We see that the system is unstable. In Table 2 are the relevant statistics for the 
nominal orbit with perturbations The performance is poorer than the nominal twodDody 
solution. The average arm length is still about 5 million km. However, the maximum 
difference between any two arm lengths is now 132,200 km and the maximum percent 
difference is now 2.634%. The maximum difference in leg rates is increased to 29.08 m/s 
and the maximum leg rate of change is 19.06 m/s. From this example, it is evident that 
we need modify the nominal orbit not only to improve the science performance, but also 
to stabilize the system for the possibility of an extended mission. 

In the next section we discuss an approach to numerically optimize the performance 
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Figure 3: Nominal Evolution with Perturbations 


of the formation using the nominal solution as an intitial guess. The purpose of the 
numerical optimization is two-fold. First, we need to improve the performance of the 
nominal orbit, and second we need to reduce the negative affects of perturbations an i 
possible, use the perturbations to our advantage. 


OPTIMIZATION APPROACH 

Here we discuss the optimization approach and details of how the cost and constraint 
functions are evaluated. The optimizer used is MATLAB’s Tmincon” functioin The 
function uses a Sequential Quadratic Progamming approach to solve the Kuhn-Tucker 


Table 2; Statistics for the Perturbed Orbit (5 Year Propagation) 


Variable Name 

2 Leg Definition 

L 

•S.OOe*^ km 

Ar 

132, 200 km 

%Ar 

2.634 % 

A?' 

29.08 m/s 

r 

19.06 m/s 
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nscessary c-'qvuitions. For a multivariablo, iion-liiioarly coiistrainod pioblom the Kuhii- 


Tucker conditions are 

v/OiO + Ea;- w;,(x*) = o ( 26 ) 

= 0 / = (27) 

A* > 0 = m.f, + 1, (28) 


where /(.x*) is the cost function evaluated at the optimum and C* is the vector of 
constraint functions evaluated at the optimum. Simply stated, the necessary conditions 
for the constrained optimal solution are that the gradient of the cost function and the 
gradients of the active constraints cancel at the solution. Furthermore, the Lagrange 
multipliers of any inactive constraints are zero at the solution point. 

VJe discussed the form of the cost functions investigated in a previous section. Here 
wc discuss some practical aspects such as what independent variables are chosen and how 
the cost functions are evaluated since analytic functions are not available. To simplify 
some of the constraints and to provide a more intuitive approach to the pioblem solution, 
we use the Keplerian elements as independent variables. Therefore, the state vector for 
the formation has 18 independent variables. For now, the epoch is assumed to be Jan 
01 2012 at zero hours. To avoid numerical problems, all angles are in radians and the 
semimajor axis is nondimensionalized on qe, the semimajor axis of the Earth s oibit 
about the sun. Lower and upper bounds are placed on the Keplerian elements to restrict 
the search space to within known feasible bounds. For example, if the semimajor axis 
is not nearly equal to 0£, then the formation will have an undesirable secular drift with 


respect to the Earth. So we require 
lower bounds applied are 

that .9a.£; 

< a < 1.1a,,. 

The complete upper and 

,9aE 

< a < 

1 . 1 a^; 

(29) 

0 

< e < 

.0105 

(30) 

0 

< i < 

.021 rad 

(31) 

0 

< CO < 

27T 

(32) 

0 

<Q < 

27T 

(33) 

0 

< M < 

27T 

(34) 




(35) 


All gradients and the Hessian matrix are approximated using finite differencing. A 
numerical approach is used to approximate the cost functions since analytic solutions are 
not available. To approximate the cost function, two years of ephemeris in one day time 
steps is used. The integrals in Eqs.(5 - 9) are approximated using a trapezoidal approach. 
Each cost function is normalized on a characteristic value to ensure the nominal solution 
is on the order of one. 
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In the next section we discuss the results of the optimization for each cost function. 
The solutions are discussed in light of LISA performance requirements and compared to 
the nominal solution. 


RESULTS 

In this section we discuss results for each of the five cost functions. A total of about 
60 converged solutions have been found. The solutions presented here are representative 
of the solutions found to date. Because there are multiple constiaints and multiple 
measures of performance, some solutions not presented here might have slightly better 
performance in one area while poorer performance in another, oi lequiie lelaxing some 
of the constraints. 

In Table 3 is statistics for the nominal orbit for a two year propagation. Ihe data 
is provided to allow a comparison between the performance of the optimal solutions and 
the nominal orbit under perturbations. 

Table 3: Statistics for the Perturbed Orbit (2 Year Propagation) 


Variable Name 

2 Leg Definition 

L 

S.OOe® km 

Ar 

61, 615 km 

%Ar 

1.232 % 

Ar 

12.072 m/s 

r 

6.918 m/s 


Results for c/i 

Recall that by optimizing according to the c/i performance measure, we are trying to 
reduce the difference in length, of two of the legs in the formation over the entire mission 
life. The orbital elements for the solution presented here are in Table 9 m Appendix 1. 
As expected the solutions all have nearly the same semimajor axis and eccentricity. Phe 
inclinations differ by about a tenth of a degree. 

In figure 4 we see the evolution of some useful parameters of the formation over two 
years. In the upper right-hand plot, it is seen that the lengths of Legl2 and Leg31 are 
nearly equal for the mission life. As expected the rate of change Legl2 and Leg31 are 
also nearly equal. It is interesting to note that two of the angles are also nearly equal m 
time. 

Table 4 shows the statistics associated with this particular optimal solution for cfj. 
The average leg length is about 5 million km. The maximum difference between the 
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X 10* Range vs. Time 



Angle vs. Time 



Range vs. Time 




Figure 4: Results for cfi 


legs of interest is 5850 km and the maximum percent change, %Ar, is .117%. This is 
an improvement of an order of magnitude over the nominal solution, for thi.s solution 
all constraints are satisfied and inactive. The formation remains about 19" behind the 
Earth for the two year mission life. 


Table 4: Statistics for Optimal Solution (c/i) 


Variable Name 

2 Leg Definition 

L 

4.982e® km 

Ar 

5850 km 

%Ar 

0.117 % 

Ar 

3.913 m/s 

r 

14.625 m/s 


Results for c/2 

By minimizing c /2 we minimize the difference in the rate-of-change of two legs of the 
formation. The orbital elements for the solution presented here are in Table 10 found 
in Appendix 1. The statistics for the solution presented here are found in Table 5. The 
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average leg length is about 5 niillion km, and is Ar is 
% and Ar and r are 2.185 m/s and 14.750 respectively, 
geometry is shown in Figure 5 


lL182km. The is 0.226 

The evolution of the formation 
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Figure 5: Results for c/2 


Table 5: Statistics for Optimal Solution (c/2) 


Variable Name 

2 Leg Definition 

L 

4.958e® km 

Ar 

11, 182 km 

%Ar 

0.226 % 

Ar 

2.185 m/s 

r 

14.750 m/s 


Results for c/3 

Table 11 contains the orbital elements for an optimal solution according to the cost 
function c/3. Recall that c/3 is developed to minimize the average rate of change of the 
legs. A figure showing the evolution of the formation geometry is in Fig. 6. The case 
shown here is for a 5 year propagation so we must compare the statistics for this case 
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shown in Table 6 with the data for the nominal orbit with perturbations for a five year 
propagation shown in Table 2. 

We see that the average leg length is about 4,93 million km. For the optimal solution, 
the maximum leg length is 61,442 km, while the maximum leg length is 132,300 for 
the nominal solution. The %Ar is 1.24% for the optimal case compared to 2.63%: for 
the nominal case. Generally speaking, for the results of the c/3 optimization, all of the 
statistics are reduced by a factor of two , as compared to the nominal solution. 


-nf Range vs. Time Range vs. Time 

X 10 ^ 






Figure 6: Results for c/3 


Table 6: Statistics for Optimal Solution (c/3) 


Variable Name 

2 Leg Definition 

L 

4.931e^ km 

At 

61442 km 

%Ar 

1.246 % 

Ar 

12.581 m/s 

r 

7.869 m/s 
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Results for c/4 and c/5 


Performance measures c/4 and c/5 are formulated to ensure all legs of the formation 
are equal in time. This is equivalent to a “breathing" equilateral triangle formation. To 
date no solutions that have acceptable performance and meet all of the design constraints 
have been found. Some solutions, which we present here, have been found by diamatically 
relaxing some of the design constraints. 

In Figure 7 we see a solution obtained using the c/5 cost function, where .several of 
the constraints have been relaxed. Specifically', we have allow^ed the leg rates to vary up 
to 200 m/s and changed the lower bound on the leg lengths to 4 million km. Although 
the figure show's tw'o years w'orth of data the optimization was only' over six months. The 
leg lengths are nearly equal for the first six months. The statistics are shown in Table 
7. The maximum difference in the leg lengths over the first 6 months is 7750 km and 
the %Ar is .160%. However, the maximum in leg rate of change is 120.25 m/s w'hich is 
much higher than allowable. Furthermore, the sy'stem is highly unstable. 



Figure 7: Results for c/5 


In Table 8 w'e see a plot of the evolution of a solution obtained using c/4. As in the 
previous example, w'e have allow’ed the leg rates to vary'’ up to 200 m/s and changed the 
low'er bound on the leg lengths to 4 million km. The average leg length for this particular 
solution is 4.2 million km. The maximum difference in leg lengths. Ar, is 34,055 km 
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Table 7: Statistics for Optimal Solution (c/5), (6 Month Optimization) 


Variable Name 

3 Leg Definition 

L 

4.843e- km 

Ar 

7750 km 

%Ar 

.160 % 

Ar 

12.16 m/s 

f 

120.25 m/s 


and %Ar is .817%. The maximum difference in leg rate of change is only 6.94 m/s, 
however, the maximum absolute leg rate is 42.02 m/s which is significantly higher than 
the acceptable value of 15 m/s. 
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Figure 8: Results for c/4 


CONCLUSIONS AND FUTURE WORK 

In this paper we presented optimal orbit configurations for the Laser Interferometer- 
Space Antenna (LISA). The performance measures used are directly related to the science 
performance. Two categories of measures were considered. The first category investigated 
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Table 8: Statistics for Optimal Solution (c/4) 


Variable Name 

3 Leg Definition 

L 

4.1C9e- km 

Ar 

34, 055 km 

<1 

0.817 % 

Ar 

6.936m/s 

r 

42.016 m/s 


functions only of two legs of the formation. The second category examined solutions that 
used all three legs of the formation in the performance metric. Several constraints were 
imposed on the .solution. The leg lengths were required to be within two percent of 
the nominal , 5 million km length. The leg ra,tes of change were constrained to be less 
than 15 m/s. Due to articulation limitations of the instrument, the interior angles of 
the forma.tion were constrained to be 60 degrees ± 1.5 degrees. The formation was also 
constrained to remain within 20 degrees of the Earth as measured from a. vertex at the 

Sun. 

Feasible optimal solutions for category one performance measures were presented. 
About one order of magnitude of improvement in the measures were achieved over the 
nominal solution. Future work in this area involves a more thorough .search of the design 
space for improved locally optimal solutions. Furthermore, relaxing the constraints may 
allow significant improvement in the perfoimance. 

It has proven difficult to find feasible optimal solutions for category two performance 
measures. Two optimal solutions were presented to illustrate possibilities if some of the 
design constraints can be relaxed. From experience gained so far, it is very difficult, if 
not Impossible, to find feasible optimal solutions to category two measures within the 
current design constraints. 
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APPENDIX 1 


Table 9: cf] Results:HeIiocentric Mean Equinox and Equator J2000 . 01 Jan 2012 00:00:00 


OE OSAl LISA2 LISA3 


a 

149563305282.0 

149563443458.0 

149563663686.0 

e 

0.009346203767 

0.009676975440 

0.0098097142755 

i 

0.9438362726 

1.000909246 

0.9196112710 

'U) 

86.37744956 

92,20962395 

91.05206384 

Q 

217.3220002 

333.4463294 

92.13655101 

u 

137.7874426 

15.41162731 

256.4720776 


Table 10: cfo Results: Heliocentric Mean Equinox and Equator J2000 , 01 Jan 2012 
00 : 00:00 - 


OE 

LISAl 

LISA2 

LISAS 

a 

152044712039.0 

152044847970.0 

152044235805.0 

e 

0.009392823346 

0.009730991627 

0.0095191410.55 

i 

0.8364363.303 

1.007377459 

1.009488404 

LJ 

90.78370457 

84.18724496 

94.71165172 

Q 

173.4799425 

299.1037480 

48.20624884 

u 

176.7385178 

.58.60738171 

297.0484844 


Table 11: c /3 Results: Heliocentric Mean Equinox and Equator J2000 , 01 Jan 2012 
00 : 00:00 - 


OE 

LISAl 

LISA2 

LISA3 

a 

149578965033.0 

149.580414155.0 

149579800648.0 

e 

0.01034342962 

0.009574892492 

0.008682661929 

i 

1.023240970 

0.9358205364 

0.8785173418 

(jj 

90.76228189 

89.11274610 

89.41129865 

n 

212.7394780 

325.7412867 

94.24688340 

u 

136.5015274 

24.79455200 

254.5655484 
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